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Abstract 



The present work is an extension of a previous study of the nucleon- 
nucleon interaction based on the chromo-dielectric soliton model. The former 
approach was static, leading to an adiabatic potential. Here we perform a 
dynamical study in the framework of the Generator Coordinate Method. In 
practice, we derive an approximate Hill- Wheeler differential equation and ob- 
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tain a local nucleon-nucleon potential as a function of a mean generator coor- 
dinate. This coordinate is related to an effective separation distance between 
the two nucleons by a Fujiwara transformation. This latter relationship is 
especially useful in studying the quark substructure of light nuclei. We inves- 
tigate the explicit contribution of the one-gluon exchange part of the six-quark 
Hamiltonian to the nucleon-nucleon potential, and we find that the dynamics 
are responsible for a significant part of the short-range N-N repulsion. 

PACS number(s): 24.85.+p, 21.30.+y, 13.75.Cs, 12.39.Ba 



2 



Typeset using REVTpjX 



I. INTRODUCTION 



In a previous investigation ffl, we studied the N-N interaction in the framework of the 
chromo-dielectric soliton model from a static point of view: we used the Born-Oppenheimer 
approximation to derive an adiabatic N-N potential, which showed a soft core repulsion due 
essentially to the color-electrostatic part of the one-gluon exchange. Previous studies of the 
N-N interaction in terms of quark degrees of freedom have pointed out the importance of 
dynamical methods (such as Generator Coordinate or Resonating Group) in the calculation 
of a realistic N-N potential. For example, in a preceding application of the non-topological 
soliton model to the N-N problem, Schuh et al. || showed that a significant part of the 
repulsion was due to dynamics; the absence of a repulsive core in some previous works was 
also interpreted as an artifact of the adiabatic approximation ||. 

The Lagrangian of the chromo-dielectric model is defined as in Ref. 

C = C q + C a + C G , (1) 

with 

C q = i> (i-fDp - m q ) V> , (2) 

C <T = ^d lt a& i a-U(a) , (3) 

£g = ~K{a)F^FIT , (4) 

where ip is the quark operator and m q the current quark mass matrix, set here to m q = 0. 
The quark Lagrangian L q is expressed in terms of the covariant derivative _D M = d il —ig s T a A° t;Ll 
and = d^Al - d v A^ + g s f abc A^A c u is the 5C/(3)-color tensor, where f abc are the SU(3) 
structure constants and T a the SU (3) generators. The quantity U(a) is the self-interaction 
of the scalar field, a, taken to be of the form: 

and the dielectric function n(a) is: 
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K (a) = 1 + 0(<r) 
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where a v is the scalar field's vacuum expectation value and 9 the usual step function. 

The quark self-energy, due to interactions with confined gluons in the dielectric medium, 
generates an effective coupling between the quarks and the scalar field: 

Cqa = -9eff(p)iH> 5 ( 7 ) 

we choose g e ff(&) to be of the form: 

9eff(<y) = go°v - 1^ • (8) 

The expression in Eq. @ is an approximation to what has been calculated in Ref. 0, and 
it is constructed to simulate spatial confinement already at the mean field level. Note that 
the coupling in Eq. (fflj breaks the chiral invariance of the Lagrangian of Eq. (fl|) . This is an 
example of dynamical symmetry breaking from which a massless Goldstone boson emerges 
naturally. 

The parameters involved in our calculation are a, b, c, go and a s = ^/47r, as discussed in 
detail in Ref. 0. By fitting the nucleon and the A masses and the proton's rms charge radius 
one remains with two free parameters, for which it is convenient to use the dimensionless 
quantities c and / = b 2 /ac. In this paper, we have chosen the set / = oo and c=10000 
taken from Table 1 of Ref. []J. Contrary to Ref. 0, the quarks here are not only coupled 
to the cx-field but also interact among themselves through one-gluon exchange (OGE). The 
OGE is treated in Abelian approximation, and it can be separated into two terms: a self- 
interaction term (in addition to g e ff(a) of Eq. (§)), which is required for color confinement 
and which contributes to the one-body part of the Hamiltonian, and a term of mutual 
interactions, which gives rise to the two-body part of the Hamiltonian. As mentioned earlier, 
in the adiabatic approximation of Ref. M, it was the color-electrostatic part of the OGE, 
which arises from the time-component of the gluonic quadrivector A c , and especially the 
corresponding self-energy diagrams, which were responsible for the soft-core repulsion. 



In this work, we incorporate the dynamics of the N-N interaction by employing the 
Generator Coordinate Method (GCM); we derive an approximate differential equation for 
the N-N wave function describing the relative motion of the two nucleons. This equation 
contains a local N-N potential and an effective, coordinate dependent mass. By means of a 
Fujiwara transformation, we then define a N-N separation length, x, from the deformation 
parameter used previously in the adiabatic approximation. This allows us to introduce a 
constant mass and to rewrite the effective potential in terms of this coordinate x. One of 
our objectives is to study the explicit role of the one-gluon exchange effects on the local 
N-N potential, included for the first time in such type of calculations. Another aim is to 
establish a connection between our effective deformation parameter and the true internucleon 
separation. The latter will enable us to apply our six-quark wave functions to studies of the 
quark substructure of light nuclei, as has been carried out already, for instance, in Ref. 0. 
The present numerical results correspond to the (TS) = (10) sector, although the formalism 
at hand can easily be extended to other isospin-spin channels. 



The GCM was introduced in the fifties by Hill and Wheeler || to describe collective mo- 
tion in nuclear systems, such as rotation, vibration or center of mass motion []7|J|. Starting 
from a many-body wave function | a ) depending on a collective coordinate a (the deforma- 
tion parameter of the six-quark system in our case), a trial wave function is constructed by 
taking a linear combination of the states | a ) with some weight function $(ct), 



II. THE GENERATOR COORDINATE METHOD 




(9) 



where $(«) is determined through the variational principle 



6E 



5 (tf|ir|tf) 



, 



(10) 



which leads to the Hill- Wheeler integral equation: 
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f(a\H -E\a) $(«') da 1 = . (11) 

This is a homogeneous Fredholm-type equation of the first kind, notoriously unstable numer- 
ically. Although some methods exist to make it stable (such as regularization || , removal of 
the zero normalization eigenmodes [|l(|, Gaussian transform fill, etc.), we prefer to solve a 
differential equation approximately equivalent to the Hill- Wheeler equation, both for numer- 
ical stability and to facilitate comparison with analyses based on the Schrodinger equation. 
In general, a is a multidimensional parameter. It is at least three-dimensional when cor- 
respondence is made to r. We here restrict the calculations to the zero-impact parameter 
case, which reduces the problem to a one-dimensional one, and leave consideration of the 
angles to a later study. 



III. THE HILL- WHEELER DIFFERENTIAL EQUATION 



To derive such a differential equation, it is more convenient to work with mean and 
relative deformation parameters, f3 and 5, defined as 

a + a' 

' (12) 
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a — a 



Expanding the weight function in a Taylor series around 5 = 0, one has: 



(V\H-E\V) = J dp J 



dd 



\(3+ 5 -\H-E\(3- 5 -) 



m - l&w + +■■■ 



(13) 



It is convenient to introduce the moments: 

H n = Jd5 (p + -\H\dp--) 5 n , (14) 
N n = JdS (p+ 6 -\dp- 6 -) 5 n . (15) 

Because ((3 + 8/2\H — E\/3 — 8/2) is an even function of 5, the odd moments are zero. 
Supposing, moreover, that (fl + 8/2\H — E\(3 — 8/2) is a sharply peaked function of 8, one 



can stop the expansion at second order in 5. Partial integration and variation by 5&* leads 
then to the Hill- Wheeler differential equation: 



1 d (,„ „„.d®\ 

__ Uh 2 -en 2 )-\ 



ld 2 H, 



dp 2 



$ = E 



\r ±d 2 N 2 



The introduction of a new function into the hermitian, 



<§(/?) = v/iVo(/3) $03) 



where 



nT 1 d 2 N 2 

allows us to transform Eq. (|16D into hermitian form: 

d 1 ci 



dp2B(p)d(3 



+ V((3) 



d>(/3) = 



where V{(3) is given by: 



N » 2M l d l } \ dfi \\jN , 



with 



Hq — H + 



ld 2 H 2 
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The term S(/3) is the effective mass: 



Nn 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



H 2 - EN 2 

The total energy E enters the definition of B; its asymptotic form at threshold is: 



(22) 



E = 2m 



N 



(23) 



where is the nucleon mass. Note that because we didn't incorporate center of mass 
corrections the asymptotic value of the potential in Eq. ( PUD is not equal to the experimental 
value of 2m n. We have indeed V(oo) = 2468 MeV when gluons are not included and V(oo) 



= 2240 MeV when they are. In practice, we could obtain a value closer to the experimental 
value by subtracting recoil corrections from the asymptotic energy: 

™^(^J-W, (24) 
but we prefer to avoid this step. This simplification does not affect our conclusions. Following 



Brink and Banerjee |L2|, we replace E in the mass term by: 



E - HM (25) 

E iVoM" (25) 

This approximation is consistent with neglecting higher order derivatives of the moments in 
the Hamiltonian. 

The moments H n (n=0, 2) and the corresponding quantities B(/3) and V(f3) have been 
calculated for three distinct cases: 

(a) H = H\ a9 + Hoge , 

(b) H = H b r + H™ 9 E , (26) 

(c) H = H\ a9 , 

where H\ a9 , Hqq E and H OGE are, respectively, the non-gluonic one-body term of the Hamil- 
tonian, the color-magnetic and the full one-gluon exchange contribution; they are given ex- 
plicitly in Ref. [0. In case (c), the one-gluon exchange was left out altogether. This is in 
the spirit of an earlier investigation where the Friedberg-Lee soliton model was applied to 
N-N scattering without considering gluonic degrees of freedom |J. In case (b), the color- 
magnetic hyperfine interaction was accounted for, and in case (a) the full color-magnetic 
and color-electrostatic OGE was included. The reason to distinguish between cases (a) and 
(b) is that in the literature it was claimed that the color-magnetic part of the OGE itself is 



responsible for the repulsive core of the N-N interaction |13[]. We shall return to this point 
at the end of Section V. 

The plot of B as a function of (3 is given in Fig. 1 for the three cases (a), (b) and 
(c). B converges towards a constant value /i, which can be calculated from considering two 
well-separated non-interacting three-quark bags: 
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H = 763.6MeV 
li = 502.3MeV 



in cases (a) and (b) 



m case c . 



(27a) 
(27b) 



We would expect \i to be equal to the reduced mass, mjv/2. The discrepancy between the 
values of /i and tun - which is especially drastic if the OGE is included, i.e., in cases (a) 
and (b) - is related to the well-known Peierls-Yoccoz disease IpjjLQl- 
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FIG. 1. The effective mass, B((3) of Eq. (|22|), as a function of the deformation parameter (3; 
the solid, dashed and dotted lines correspond, respectively, to the cases (a), (b) and (c) introduced 



in Eq. (26). The asymptotic values of Eqs. (27a) and (27b) are indicated by the arrows. 



IV. THE FUJIWARA TRANSFORMATION 

The dependence of the effective mass on f3 prevents us from directly interpreting the 
potential in Eq. (p0| ) as an ordinary N-N potential. Moreover, (3 doesn't correspond to 
the true N-N separation distance (except for large positive j3 when the two nucleons are 
well separated). Therefore, we wish to transform Eq. (|TjJ) into a Schrodinger-like equation 
with a constant, coordinate independent mass term. For this purpose, one can use a Fuji- 
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wara transformation [|14||15|| , which relates the generator coordinate {3 to an effective N-N 
separation length: 



x(f3) 



d/3' + /3. 



If one now redefines the weight function in Eq. (|19|) as 



Eq. (|19|) transforms into the familiar form 

1 d 2 



2/x dx 2 

with V given by Eq. (|20|) and 



v + v F 



ijj(x) = Ei])(x) , 



7 (dB\ 2 _ 1 rf 2 5 
32fi 3 ~ 8B 2 d(5 2 
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FIG. 2. The Fujiwara coordinate, x{0) of Eq. (]28|), as a function of the deformation parameter 
j3. The long-dashed line shows the asymptotic solution, x{(5) = (3, and the remaining labeling is 
the same as in Fig. 1. 

Figure 2 displays the explicit relationship between x and /3, as obtained from Eq. ([Z8"D. 
As expected, the deformation parameter (3 converges asymptotically towards the effective 
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internucleon separation x. The correspondence (3 <-> x should be very useful in discussions 
of the quark substructure of nuclei or nuclear matter using Schrodinger-based many-nucleon 
calculations and employing our six-quark wave functions. 



V. RESULTS FOR THE EFFECTIVE N-N POTENTIAL 



We now wish to present detailed results for: 



VioM = V(/3) + V F {(3) - V{oo) 



(32) 



where V((3) and Vf(P) are given in Eqs. (|20|) and (|3T|) . The value of V(oo) corresponds 
to the asymptotic value of Hq/Nq calculated from two well-separated non-interacting three- 
quark bags, and it is given in Section III. This asymptotic value is the same in cases (a) and 
(b) because the color-electrostatic mutual and self-energy terms cancel exactly due to color 
neutrality when the two nucleons are well separated. 

It is convenient to rewrite Vi oc (P) in the following form: 



V loc (P) = V (P)+V 1 (P) 



(33) 



with: 



v {p) 



N Q 



i 

4B 



d 2 In N 1 

d(3 2 + 2 



'd\nN n 



dp 



d\nN d\nB 

~df3 W 



1 

8B 



3 f dlnB ' 
4 



d 2 In B 

dp 2 



(34) 



(35) 



In order to calculate these derivatives, InB and IniVo were fitted to polynomials. The two 
contributions Vi(ft) (i=0, 1) to Vi oc (P) are plotted in Fig. 3 as functions of the deformation 
parameter P for the three cases outlined previously. Fig. 4 shows Vi oc as a function of P and 
of the Fujiwara coordinate x, respectively. Note that Eq. ( |35| ) was obtained from Eq. ( |20|) 
by replacing H 2 — EN 2 with —N /B, as indicated in Eq. (|22|). 
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FIG. 3. The two contributions to the local potential, Vq((5) of Eq. ( |3jJ ) and Vi(f3) of Eq. 
as functions of the mean generator coordinate (3. The solid, dashed and dotted lines correspond to 
the cases (a), (b) and (c) introduced in Eq. (p6|). 
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FIG. 4. The non-adiabatic, local potential, Vi oc of Eq. (p3|), as a function of the deformation 
parameter (3 and the Fujiwara coordinate x, respectively. The labeling is the same as in Fig. 3. 
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The shape of Vq(/3) is quite similar to the adiabatic potentials displayed in Fig. 10 of Ref. 
|Ij], both for the "full OGE" and "no-OGE" cases. This tends to confirm our assumption 
that the matrix elements (/3 + 5/2\H — E\/3 — 5/2) are rather sharply peaked around 5 = 0. 
The term V\ (/?) corresponds to the contribution of non-adiabaticity. It grows important only 
for (3 < — 2 fm, and yields in all cases a repulsion due to the dynamics. This is according to 
our expectation and in agreement with Ref. ||. Note that in cases (b) and especially (c), 
we also obtain an intermediate range attraction in V/ oc . The fact that our N-N potential 
extends to negative x should not be taken too literally. It simply reflects inadequacies in the 
relationship between the deformation parameter j3 and the N-N separation length x, which 
are connected to the Peierls-Yoccoz disease mentioned earlier. 

We recall that one of the main objectives of this and our previous study |J was to 
incorporate explicitly one-gluon exchange effects, in contrast to Ref. [Q] where they were 
neglected. Comparing, for instance, cases (a) and (c), one can see that the OGE reinforces 
the repulsive core considerably. The existence of a repulsive core in all three cases makes us 
to attribute it to dynamics rather than to the color-magnetic interaction (case (b)), as was 
inferred in Ref. [|Hf|. 



VI. CONCLUSIONS 

In this investigation, we found that the dynamics are manifestly responsible for the 
hard-core repulsion of the short-range part of the N-N interaction, and we observed that we 
could obtain both short-range repulsion and some intermediate range attraction if the entire 
one-gluon exchange or at least its electrostatic part were neglected. 

In the results containing the full OGE effects the lack of attraction is due to the omission 
of explicit meson exchanges. Then, to reproduce the experimental phase shifts or other two- 
body data one necessitates to attach a local OBE potential beyond a certain internuclear 



distance To obtain this potential in the framework of our model we could consider 



extending our calculations by either including quantum surface fluctuations and/or intro- 
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during configurations of the form q 7 q in addition to the q 6 states. This would be a rather 
cumbersome procedure within the present model. The most convenient would be to either 
allow mesonic degrees of freedom and to consider, e.g., an explicit pion exchange between 



the individual quarks (T7j or to simply choose a phenomenological potential. 

Another important result of this work is the evaluation of the relationship between the 
deformation parameter (3 and the effective N-N separation length x through the Fujiwara 
transformation. This correspondence is very useful for applications of our model to the 
description of phenomena involving the quark substructure of light nuclei. It furthermore 
allows us to relate many-body correlation functions or N-N wave functions given in the 
literature to the GCM formalism presented here. 

An attractive way to confirm our results would be to solve directly the Hill- Wheeler 
integral equation in order to obtain phase shifts. Projection on good angular momentum 
states should also improve our calculations. 
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